Heat transport in quantum spin chains: the relevance of integrability 



Jinshan Wu and Mona Berciu 
Department of Physics and Astronomy, University of British Columbia, Vancouver, BC, Canada, V6T 1Z1 

(Dated: March 9, 2010) 

We investigate heat transport in various quantum spin chains, using the projector operator tech- 
nique. We find that anomalous heat transport is linked not to the integrability of the Hamiltonian, 
but to whether it can be mapped to a model of non-interacting fermions. Our results also suggest 
how seemingly anomalous transport may occur at low temperatures in a much wider class of models. 
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Heat transport in quantum spin chains, in particular 
when is normal (diffusive) transport observed, is still not 
understood despite considerable effort^— For example, 
it was conjectured that integrability leads to anomalous 
(ballistic) transport^ but it was also argued that an inte- 
grable gaped XXZ chain has normal conductivity^ Oth- 
ers have argued that only the spin conductivity is nor- 
mal in this case, while the thermal conductivity is still 
anomalous A consensus on what are the necessary cri- 
teria for normal conductivity is still missing. 

Most of the above worker— studied infinite and/or pe- 
riodic chains, and used the Kubo formul a 12 ' 13 where fi- 
nite/zero Drude weight signals anomalous/normal trans- 
port. For integrable systems, the Kubo formula always 
predicts anomalous heat transport. In fact, full integra- 
bility is not even necessary, all that is needed is com- 
mutation of the heat current operator with the total 
Hamiltonian. 4 Anomalous heat transport observed ex- 
perimentally in systems described by integrable models, 
such as (Sr,Ca)i 4 Cu240 4 i, Sr 2 Cu0 3 and CuGe0 3 r^ 
seems to validate this result, although Ref. Qj] finds nor- 
mal transport in Sr 2 Cu03 at high temperatures. 

The proof of the Kubo formula requires dealing with 
currents between the chain's ends and the thermal baths 
it is connected to.— For infinite systems, one may argue 
that such currents can be ignored, as they are a bound- 
ary effect. However, the terms describing the coupling to 
the baths lead to a non-vanishing commutator between 
the heat current operator and the total Hamiltonian, in- 
validating the main argument for anomalous transport. 
In other words, "integrability" of the chain connected to 
baths may be lost even if the isolated chain is integrable. 
This conclusion is supported by recent proofs of Kubo- 
type formulae for finite systems, based on phenomeno- 
logical approaches such as the Fokker-Planck equation. 19 
These Kubo formulae have similar structure to the orig- 
inal one, however the dynamics is not defined only by 
Hamiltonian of the chain but also includes random vari- 
ables mimicking the effects of coupling to the baths. 

Here we investigate finite spin chains coupled to ther- 
mal baths, using the projector technique^ - — There are 
many other similar studie o 11 ' 21 " — for various spin Hamil- 
tonians, most of which would be integrable for a periodic, 
isolated chain (hereafter we call such models integrable). 
Results range from normal to anomalous transport, and 
there is no agreement on whether integrability is corre- 



lated or not with anomalous transport. 

We propose a resolution for this question in this Rapid 
Communication. We find that integrability is not a suf- 
ficient condition for anomalous heat transport. We find 
anomalous transport at all temperatures only in models 
which can be mapped onto homogeneous non-interacting 
fermionic models. All other models we investigated ex- 
hibit normal heat transport, whether they are integrable 
or not (however, as discussed below, at low temperatures 
their heat transport may become anomalous in certain 
conditions). We therefore conjecture that the existence 
of such a mapping is the criterion determining anomalous 
transport, at least for finite-size systems. 

We begin by briefly describing our calculation method, 
which is a direct generalization of the projection operator 
technique used to study the evolution towards equilib- 
rium of a system coupled to a single bath^ The iV-site 
chain of spins- h is described by the Hamiltonian: 

N-l N 
i=l i=l 

while the heat baths are collections of bosonic modes: 

HB = y^^k,gb{. a bk.a 

where a = R/L indexes the right/left-side baths and we 
set fi = l,ks = 1, and the lattice constant a = 1. The 
system-baths coupling is taken as: 

k,a 

where ir, — 1 and iji = N, i.e. the left (right) thermal 
bath is only coupled to the first (last) spin and can induce 
its spin-flipping. This is because we choose B ■ e y = 
while \B\ is finite, meaning that spins primarily lie in the 
xOz plane so that s y acts as a spin-flip operator. 

The evolution of the total system is described by the 
Liouville-von Neumann equation for the total density ma- 
trix pt- If we are interested only in the properties of 
the central system, it is convenient to find an equation 
of motion for the reduced density matrix ps — tr B (pt) 
and solve it directly. This is achieved by using the pro- 
jection operators, treating the system-bath coupling per- 
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turbationally to second order, and also by using a Marko- 
vian approximation! 20 ! 21 These approximations are rea- 
sonable: the system-baths coupling must be rather weak 
so that the properties of the chain are determined by its 
specific Hamiltonian, not by this coupling. The Marko- 
vian approximation is also justified, since we are inter- 
ested only in the steady-state limit t — > oo. 
The resulting equation of motion for ps (t) is: 

dpsjt) 
dt 



= -i{H s ,Ps(t)}-\ 2 J2 ([ s i a >™<xPs(t)]+h.c.) 



a=L,R 

(1) 

where rh a = sf^ ■ S Q . Here, (•) refers to the element-wise 
product of two matrices, (n\a ■ b\m) — (n\a\m) (n\b\m) . 
The bath matrices £l,.r are defined in terms of the eigen- 
statcs of the system's Hamiltonian Hs\n) = E n \n) as: 



t a = n*^2 \m)(n\ 6 (0 TOn )n Q (fi mn ) D a (fi m „) \ V h 



(a) |2 



(<*) |2 



-Q (Cl nm ) (1 + n a (fi nm )) D a (n nm ) \V£*> | 



and L„ is defined 



where Q, r , 

transition. Furthermore, Q(x) is the Heaviside function, 



= n r 



a bath mode resonant with this 



n a (Q.) 



-1 



is the Bose-Einstein equilibrium 



distribution for the bosonic modes of energy Q at the 
bath temperature T a = l//3 a , and D a (£l) is the bath's 

density of states. The product D a (Ct mn ) \v£ | 2 is the 
bath's spectral density function. For simplicity, we take 
it to be a constant independent of m and n. 

The sam e eq uation of motion for ps{t) was also derived 
in Refs. [2l| - |23l where its steady-state solution was found 
via Runge-Kutta integration or by solving an eigenvalue 
problem. The latter comes about because the steady- 
state poo is given by = Cpoo, where L is the linear 
operator on the RHS of Eq. ([T]) , so in matrix terms poo 
is the eigenvector corresponding to the zero eigenvalue 
of C. Using the normalization trpoo = 1, this eigenvalue 
problem can be replaced with solving a linear system 
of coupled equations, which makes it more efficient and 
allows us to analyze somewhat larger systems. 

We rewrite Hs = YliLl 1 fh,i+l+Y^=i where h l4+1 
is the exchange between nearest-neighbor spins and hi 
is the on-site coupling to the magnetic field. We can 
then define a local site Hamiltonian h 
hi + hhi,i+i (with ho t 



(S) 



- 



Hamiltonian h 



(B) 



h n,n+i = 0) and a local bond 
- hi t i+i + \hi+i such that %s = 

J2iLi ~ SiLi 1 ^i^ - The local site Hamiltonians 
can be used to derive the heat current operator from the 



2 hi 



continuity equation ji 



Ji-l- 



= Vj = - 



dt 



H s ,h 



(B) 



h (B) h (B) 



for 



This results in ji^i+i = i 

■ , N — 2. As expected, in the steady state we 
= J to be independent of i. 



1 = h 
find (j. 

Knowledge of the steady state heat current J, as such, 
is not enough to decide whether the transport is normal 



or not. Consider an analogy with charge transport in a 
metal connected to two biased leads. What shows if the 
transport is anomalous is the profile of the electric poten- 
tial, not the value of the electric current. In anomalous 
transport (clean, non-interacting metal) all the voltage 
drop occurs at the ends of the sample, near the contacts. 
Away from these contact regions, electrons move ballisti- 
cally and the electric potential is constant, implying zero 
intrinsic resistance. For a dirty metal, scattering takes 
place everywhere inside the sample and the electric po- 
tential decreases monotonically in between the contact 
regions, i.e. the sample has finite intrinsic resistivity. 

In principle, the scaling of the current with the sam- 
ple size, for a fixed effective bias, also reveals the type 
of transport: for anomalous transport, the current is in- 
dependent of the sample size once its length exceeds the 
sum of the two contact regions, while for normal trans- 
port it decreases like inverse length. The problem is that 
one needs to fix the effective bias, i.e. the difference be- 
tween the applied bias and that in the contact regions. 
Furthermore, since we can only study relatively short 
chains, the results of such scaling may be questionable. 

It is therefore desirable to use the equivalent of the 
electric potential for heat transport and to calculate its 
profile along in order to determine the type of transport. 
This, of course, is the "local temperature", which is a dif- 
ficult quantity to define. One consistency condition for 
any definition is that if Tl = Tr = T, i.e. the system 
is in thermal equilibrium at T, then all local tempera- 
tures should equal T. We define local site temperatures 
Ti which fulfill this condition in the following way. Since 
we know all eigenstates of Hs, it is straightforward to 
calculate its equilibrium density matrix at a given T, 



Ps 



eq,T _ 1 



Z 



\n){n\, where Z = Y. n e ~ PEn - Lct 



then (h[ S) ) eq _ T = tr[p e ^ T hf\ We define T f to be the 

(S) ( s) 

solution of the equation: (h\ ) e q,Ti — tr[pooh\ ]. In 

other words, the steady-state value of the energy at that 

site equals the energy the site would have if the whole 

system was in equilibrium at Tj. Of course, we can also 

( B) 

use other "local" operators such as h\ to calculate a 
local bond temperature T i+ i. We find that when these 
definitions are meaningful, the results are in very good 
agreement no matter what "local" operator is used. 

This type of definition of Ti is meaningful only if a large 
magnetic field B is applied. For small B, the expecta- 
tion values (h\ ) e q,T are very weakly T-dependent, so 
that tiny numerical errors in the steady-state value can 
lead to huge variations in T"j. Addition of a large B is 
needed to obtain (h\ ) e q,T which varies fast enough with 
T for values of interest so that a meaningful Tj can be ex- 
tracted. Since we could not find a meaningful definition 
for Ti when B — > 0, we cannot investigate such cases. 
Note, however, that most integrable models remain inte- 
grable under addition of an external field B = Be z . 

In all of our calculations, we take B z = 1 and the 
exchange J ~ 0.1. Temperatures T L / R — T(l ± 6/2) 
should not be so large that the steady state is insensitive 
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FIG. 1: (color online) Plot of (a) local temperature profile 
for chains with N = 9, and (b) thermal current J vs. chain 
length N. In all cases Tl = 2.4, Tr = 1.6. Other non-zero 
parameters are J x — J y = 0.1, B z — 1 (XY); J x = J y = 
\j z = 0.1, B, = 1 (XXZ); J x = \ J y = \.J Z = 0.1,5, = 1 
(XYZ); J z = 0.2, B x = 1 (Ising*) and J z = 0.2, B x =B Z = 1 
(Ising xz ). The XY and the Ising^ chains show flat T profiles 
and currents independent of the chain's length, i.e. anoma- 
lous transport. All other models have normal transport. 



to the model or so small that only the ground-state is 
activated. Reasonable choices lie between min( J x , J y , J z ) 
and NB, which are roughly the smallest, respectively the 
largest energy scales for an iV-site spin chain. 

In Fig. QJa) we show typical results for local tem- 
perature profiles Tj, T i+ i. We apply a large bias S = 
(Tl — Tr)/T = 0.4 for clarity, but we find similar results 
for smaller 5 (see below). For these values, the "contact 
regions" include about two spins on either end. The Ti 
profile of the rest of the chain is consistent with anoma- 
lous transport (flat Ti profile) for the XY chain and 
shows normal transport (roughly linear Ti profile) in all 
the other non-Ising, J z ^ cases. We find similar results 
(not shown) for ferromagnetic couplings. All these are 
integrable models. The XY model is special because it 
can be mapped to non-interacting spinless fermions with 
the Jordan- Wigner transformation^ A finite J z leads to 
nearest-neighbor interactions between fermions. Eigen- 
modes for J z ^ can be found using Bethc's ansatz, but 
they cannot be mapped to non-interacting fermions. 

Another model that maps to non-interacting spinless 
fermions is the Ising model in a transverse field B x *& For 
this model we again find anomalous transport, as shown 
in Fig. HJa). If we add a B z field, the model becomes 
non-integrable^ and we recover normal transport. The 
scaling of J vs. N, shown in Fig. Hfb), supports these 
conclusions, although a quantitative analysis is difficult 
because of the contact regions's contributions. 

We found this generic behavior for a wide range of 
parameters. When A G [0.03,0.2], T G [0.3,30.0] and 
5 > 0.01, the XYZ model has normal conductivity when 
J z G [0.03, 0.5] and anomalous conductivity when J z = 0. 
When < J z < 0.03 or J z G [0.5, 1] the local temper- 
ature still decreases monotonically but not linearly, and 
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FIG. 2: (a) Slope S vs. system size N, and (b) normalized 
temperature Tj ie H = (Ti — Tn-i)/(T2 — Tn-i) for various 
systems sizes N. Parameters are J x — J y = 0.1, J z = 0.2, 
B z = 1.0, T = 2.0, S = 0.2, A = 0.1. See text for more details. 



the current decreases more slowly than 1/N. Since we 
cannot study much longer chains we cannot easily distin- 
guish here between normal vs. anomalous transport. For 
J z > 1 the system becomes Ising-like and the transport 
is anomalous, as expected. 

We would like to better gauge how things change with 
N, even with our limited range. For this, we consider how 
the effective temperature bias on the chain, T 2 ~ Tn-U 
or its effective slope S = ^T^" 1 , depend on the system 
size N. For normal transport we expect S cx (N—^)^ 1 (f3 
accounts for the contact regions). If S cx (N — (3)~ a with 
a > 1, then for longer chains the temperature profile 
tends to be flatter than normal so we will use a > 1 
as a signature of anomalous transport (of course, if a 
plateau starts to emerge near the center of chain that 
also indicates anomalous transport). A second gauge of 
the size dependence comes from looking at how the shape 
of the normalized temperature profile changes with N. 

Figures^a) and (b) show such analysis. The left panel 
shows fits for S (solid lines) for XXZ and XYZ models 
(symbols). Best fits give a < 1, consistent with nor- 
mal transport. Similarly, the right panel which plots the 
normalized temperature profiles for different values of N 
shows no change with increasing N, and no evidence that 
a plateau may evolve. Based on this limited evidence, we 
conclude that these models, although integrable, do ex- 
hibit normal transport. 

In summary, the first conclusion we draw from these 
results is that integrability is not sufficient to guar- 
antee anomalous transport: several integrable mod- 
els show normal heat transport, in agreement with 
other studies i^ ' 23 ' 24 i 26 The second conclusion is that only 
models that map onto Hamiltonians of non-interacting 
fermions exhibit anomalous heat transport. This is a rea- 
sonable sufficient condition, since once inside the sample 
(past the contact regions) such fermions propagate bal- 
listically. However, we cannot, at this stage, demonstrate 
that this is a necessary condition as well. We therefore 
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can only conjecture that this is the criterion determining 
whether the heat transport is anomalous. 

In this context, it is important to emphasize again the 
essential role played by the connection to the baths. In 
its absence, an isolated integrable model is described by 
Bethe ansatz type wavefunctions. Diffusion is impossible 
since the conservation of momentum and energy guaran- 
tees that, upon scattering, pairs of fermions either keep 
or interchange their momenta. For a system connected 
to baths, however, fermions are continuously exchanged 
with the baths, and the survival of a Bethe ansatz type 
of wavefunction becomes impossible. In fact, even the 
total momentum is no longer a good quantum number. 
We believe that this explains why normal transport in 
systems mapping to interacting fermions is plausible. 

Normal transport is also possible for non-interacting 
fermions, if they are subject to elastic scattering on dis- 
order. This can be realized, for example, by adding to 
the XY model a random field B z at various sites. We 
have verified (not shown) that a local drop in the local 
temperature indeed arises near sites with such disorder, 
leading to normal conductance in "dirty" samples. 

On the other hand, anomalous transport can also oc- 
cur in models which map to homogeneous interacting 
fermions if the bath temperatures are very low. Specifi- 
cally, consider the XXZ models. Because of the large B z 
we use, the ground-state of the isolated chain is ferromag- 
netic with all spins up. The first manifold of low-energy 
eigenstates have one spin flipped (single magnon states), 
followed by states with two spins flipped (two magnon 
states), etc. The separation between these manifolds is 
roughly B z , although because of the exchange terms each 
manifold has a fairly considerable spread in energies and 
usually overlaps partially with other manifolds. 



If both Tl , Tji <C B z , only single-magnon states par- 
ticipate in the transport. We can then study numerically 
very long chains by assuming that the steady-state ma- 
trix elements p nm vanish for all other eigenstates (S z ,tot 
is a good quantum number for these models). In this case 
we find anomalous transport for all models, whether inte- 
grable or not. This is reasonable, since the lone magnon 
(fermion) injected on the chain has nothing else to inter- 
act with, so it must propagate ballistically. 

We can repeat this restricted calculation by includ- 
ing the two-magnon, three-magnon, etc. manifolds in 
the computation. As expected, the results agree at low 
Tl,Tr, but differences appear for higher Tl,Tr, when 
these higher-energy manifolds become thermally acti- 
vated. In such cases, the transport becomes normal for 
the models mapping to interacting fermions as soon as 
the probability to be in the two (or more) magnon sector 
becomes finite. In other words, as soon as multiple exci- 
tations (fermions) are simultaneously on the chain, and 
inelastic scattering between them becomes possible. 

These results may explain the heat transport observed 
experimentally in compounds such as Sr2Cu03^ where 
at low temperature anomalous transport was found while 
at high temperature normal transport was reported. 

In conclusion, we propose a new conjecture for what 
determines the appearance of anomalous heat transport 
at all temperatures in spin chains. Unlike previous sug- 
gestions linking it to the integrability of the Hamilto- 
nian or existence of gaps, we propose that the criterion 
is the mapping of the Hamiltonian onto a model of non- 
interacting fermions without any disorder. 
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